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ABSTRACT:  We  present  the  first  temperature-dependent  rate  constants  and 
branching  ratios  for  the  reactions  of  Fe+  and  FeO+  with  C2H2,  C2H4,  and  C2H6  from 
170  to  700  K.  Fe+  is  observed  to  react  only  by  association  with  the  three  hydrocarbons, 
with  temperature  dependencies  of  T~2  to  T-3.  FeO+  reacts  with  C2H2  and  C2H4  at  the 
collision  rate  over  the  temperature  range,  and  their  respective  product  branchings  show 
similar  temperature  dependences.  In  contrast,  the  reaction  with  ethane  is  collisional  at 
170  K  but  varies  as  T-0'5,  while  the  product  branching  remains  essentially  flat  with 
temperature.  These  variations  in  reactivity  are  discussed  in  terms  of  the  published 
reactive  potential  surfaces.  The  effectiveness  of  Fe+  as  an  oxygen-transfer  catalyst  toward 
the  three  hydrocarbons  is  also  discussed. 


■  INTRODUCTION 

The  catalytic  role  of  transition-metal  ions  in  the  activation  of 
small  molecules  has  garnered  significant  attention  over  the 
years.1-17  In  this  light,  gas-phase  studies  of  ion— molecule 
reactions  have  led  to  many  insights  into  the  underlying 
mechanisms  of  numerous  reactions  involving  transition  metals, 
such  as  the  concept  of  two-state  reactivity.  ’  ’  '  While  direct 

application  of  gas-phase  studies  to  industrial  processes  is 
limited,  the  fundamental  insight  into  the  underlying  reactions 
allows  for  a  more  rational  and  efficient  development  of  catalysts 
applicable  to  the  bulk  phase.  A  prime  example  of  this  is  the 
large-scale  production  of  HCN  from  methane  and  ammonia, 
the  understanding  of  which  has  greatly  benefited  from  j;as- 
phase  studies  of  Pt+  with  both  methane  and  ammonia. 

Also  important  is  the  ability  of  gas-phase  studies  to  lead  to  the 
development  of  efficient  quantum  chemical  calculation 
methods  by  offering  benchmarks  for  testing  and  refinement. 
Due  to  the  significant  computational  scaling  with  respect  to 
treatment  of  electron  correlation  for  systems  of  this  size,  it  is  an 
important  goal  to  find  an  acceptable  compromise  between 
accuracy  and  computational  effort.  The  most  commonly 
employed  technique  for  transition-metal-containing  systems, 
density  functional  theory  (DFT),  has  performed  with  mixed 
results.  Wide  variations  in  calculated  values  are  found  to  be 
highly  dependent  on  thr  functional  and  basis  set  employed, 
with  little  correlation  to  the  system  studied.”7  30 

We  present  here,  for  the  first  time,  temperature-dependent 
kinetic  data  for  the  reactions  of  Fe+  and  FeO+  with  acetylene, 
ethylene,  and  ethane.  Catalytic  activation  of  C— C  and  C— H 
bonds  is  of  considerable  interest,  and  numerous  studies  have 
considered  the  mechanistic  role  of  transition  metals  and  their 


oxides  in  these  reactions,2’12’31-38  especially  iron.21’39-48  Here, 
we  present  the  temperature  dependence  of  the  reaction  rate 
constants  and  the  product  distributions  in  an  effort  to  not  only 
gain  fundamental  insight  into  the  reaction  mechanisms  but 
additionally  to  allow  for  the  assessment  and  refinement  of  the 
calculated  potential  energy  surface  for  the  reaction  profiles. 

■  EXPERIMENTAL  METHODS 

All  measurements  were  performed  on  the  Air  Force  Research 
Laboratory’s  variable  temperature  selected  ion  flow  tube  (VT- 
SIFT)  instrument,  which  has  been  described  in  detail 
elsewhere.49  Briefly,  Fe+  ions  are  created  using  an  electron 
impact  source  in  the  presence  of  a  10%  mixture  of  Fe(CO)5  in 
He,  with  N20  additionally  supplied  in  the  source  for  the 
formation  of  FeO+.  Ions  are  extracted  and  injected  into  a 
quadrupole  mass  filter  where  the  desired  ion,  either  Fe+  or 
FeO+,  is  mass  selected.  The  ions  are  focused  before 
introduction  into  a  laminar  flow  tube  via  a  Venturi  inlet, 
where  ~104— 105  collisions  with  a  He  buffer  gas  act  to 
thermalize  the  ions  and  carry  them  downstream.  Previous  work 
in  our  lab  has  shown  that  some  excited  electronic  states  of  Fe+ 
are  inefficiently  quenched50  and  may  be  present  in  the  flow 
tube.  FeO+,  however,  is  expected  to  be  fully  quenched  due  to  its 
extra  degrees  of  freedom.  In  either  case,  linear  decays  in  parent 
ion  signal  on  a  semilogarithmic  scale  were  observed  over  several 
orders  of  magnitude,  indicative  of  a  single  rate  constant  (see  the 
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Supporting  Information).  The  temperature  of  the  flow  tube  is 
variable  over  a  large  range  either  by  resistive  heating  devices 
(300—700  K),  pulsed  liquid  nitrogen  (120—220  K),  or 
recirculating  methanol  chillers  (220—300  K).  Operating 
pressures  of  0.4  Torr  are  maintained  in  the  flow  tube 
throughout  all  experiments,  other  than  for  verification  of 
ternary  reaction  channels.  The  neutral  reagent  (C2H2,  C2H4,  or 
C2H6)  is  added  59  cm  upstream  of  the  end  of  the  flow  tube, 
with  typical  reaction  times  on  the  order  of  4  ms,  dependent  on 
helium  buffer  flow  (varied  from  10  to  13  std.  L  min-1)  and 
temperature. 

After  traveling  the  length  of  the  flow  tube,  the  core  of  the 
flow  is  sampled  through  a  truncated  nose  cone  with  a  2  mm 
aperture.  The  remainder  of  the  flow  is  pumped  away  by  a  Roots 
pump  through  a  throttled  gate  valve  that  acts  to  maintain  the 
desired  pressure  within  the  flow  tube.  After  the  nose  cone,  the 
primary  ions  and  product  ions  are  guided  by  a  lens  stack  to  a 
quadrupole  mass  filter  for  analysis  and  are  subsequently 
detected  using  an  electron  multiplier.  Rate  constants  are 
derived  by  monitoring  the  decay  of  the  primary  ion  as  a 
function  of  the  neutral  reagent  flow.  Product  branching  is 
determined  by  monitoring  the  relative  product  formation  as  a 
function  of  neutral  reagent  and  extrapolating  to  zero  flow. 
Measurements  were  made  from  approximately  120  to  700  K, 
and  temperature  dependences  of  rate  constants  and  product 
branching  were  determined.  Errors  in  the  rate  constants  are 
estimated  to  be  ±25%  absolute  and  ±15%  relative  to  each 
other.51  Errors  in  the  product  branching  are  estimated  to  be 
±10%. 

■  RESULTS  AND  DISCUSSION 

1 .  Fe+  +  C2Hx.  The  effective  tertiary  rate  constants  as  a 
function  of  temperature  for  the  reactions  of  Fe+  with  C2H2, 
C2H4,  and  C2H6  in  0.4  Torr  of  He  are  shown  in  Figure  1  and 
listed  in  Table  1.  In  all  cases,  the  only  primary  product  channel 
observed  was  association  (e.g.,  the  C2H2  reaction  yielding 
FeC2H2+).  The  association  complexes  were  observed  to  further 
react  by  association,  as  has  been  previously  reported.44  The 
reaction  of  Fe+  with  ethane  is  known  to  have  an  exothermic 
channel  producing  FeC2H4+  +  H2,39  but  this  channel  was  not 
observed  in  this  work  or  in  previous  SIFT  work  by  Bohme  et 
al.44  This  is  unsurprising  as  the  previously  reported  rate 
constant  (~l/500  the  collision  rate  constant)  is  at  the  lower 
end  of  our  detection  limit,  and  the  pressures  employed  both  in 
this  work  and  that  by  Bohme  et  al.  favor  the  association 
complex  much  more  strongly  than  did  the  low-pressure  beam 
experiments  where  that  channel  was  observed. 

Figure  2  shows  the  number  density  dependences  of  the  rate 
constants  at  300  K  over  the  range  accessible  in  the  VT-SIFT. 
Least  squares  lines  forced  through  zero  are  shown.  The 
reactions  of  both  acetylene  and  ethane  show  that  that  is  correct 
within  error,  indicating  that  they  are  in  the  low-pressure 
limit,52’53  and  the  termolecular  rate  constant  is  obtained  from 
the  slope.  In  contrast,  the  faster  reaction  with  ethylene  showed 
that  forcing  a  zero  intercept  is  a  poor  representation  of  the  data, 
likely  indicating  that  the  reaction  is  intermediate  between  the 
low-  (zero  intercept)  and  high-pressure  limits  (flat  with  P). 
Because  the  observed  C2H4  rate  constants  approach  ~10%  of 
the  collision  rate  constant  (compared  to  just  1—2%  for  the 
other  reactions),  perhaps  this  is  not  surprising.  The  alternative 
explanation  for  the  nonzero  intercept  is  that  there  is  either  a 
bimolecular  or  radiative  association  component  to  the  reaction. 
Because  no  bimolecular  products  were  observed,  the  former 
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Figure  1.  Tertiary  rate  constants  for  the  reactions  of  Fe+  with  C2H2J 
C2H4,  and  C2H6  in  0.4  T orr  of  He.  Pressure  dependences  show  that 
the  C2H2  and  C2H6  values  are  true  tertiary  rate  constants  and  that  of 
C2H4  is  an  effective  one  in  the  falloff  regime.  Relative  15%  error  bars 
are  shown. 


Table  1.  Effective  Tertiary  Rate  Constants  as  a  Function  of 
Temperature,  k(T),  for  the  Reactions  of  Fe+  with  C2H2, 
C2H4,  and  C2H6  in  0.4  Torr  of  He" 

reaction  fc(T)  cm6  S_1 

Fe+  +  C2H2  Fe(C2H2)+  1.32  X  10-27(T/300)-28±02 

Fe+  +  C2H4  Fe(C2H4)+  7.25  X  l(r27(T/300)-2-1±0-2 

Fe+  +  C2H6  -s-  Fe(C2H6)+  1.30  X  10-27(T/300)-29±0-2 

“The  C2H2  and  C2H6  values  are  true  termolecular  values,  while  the 
C2H4  values  are  effective  rate  constants  (see  the  text).  Errors  in  the 
rate  constants  are  estimated  to  be  ±25%  absolute  and  ±15%  relative  to 
each  other.51 


can  be  ruled  out.  The  observed  intercept  would  require  an 
extremely  large  radiative  rate,  which  is  unlikely.  Excited  states  of 
Fe+  can  also  be  ruled  out  as  they  would  result  in  curvature  in 
the  decay,  which  is  not  observed  over  1—3  orders  of  magnitude 
(see  the  Supporting  Information)  or  different  products.  Thus, 
the  reported  rate  constants  for  C2H2  and  C2H6  are  true 
termolecular  values,  and  those  for  C2H4  are  effective  values. 
Unfortunately,  we  were  unable  to  explore  these  reactions  over  a 
large  enough  pressure  range  to  allow  for  precise  modeling  of 
the  falloff  curves. 

Stationary  points  on  the  potential  surfaces  for  all  three 
reactions  were  calculated  using  DFT.  While  energy  ordering  of 
the  various  spin  states  is  a  well-known  issue  with  density 
functional  calculations  on  transition-metal  systems,  the 
calculations  here  employing  B3LYP/TZVP54'5'  using  the 
Gaussian  09  program  suite5  yield  a  6D— 4F  energy  splitting 
for  Fe+  of  0.29  eV,  in  reasonable  agreement  with  the 
experimental  spin— orbit  average  of  0.247  eV.57,58  In  agreement 
with  previous  calculations,59-62  the  ground  state  of  each 
association  complex  was  found  to  be  quartet,  as  opposed  to  the 
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Fe+  +  CH  ->  Fe(C  H  )+ 
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Figure  2.  Effective  bimolecular  rate  constant  for  the  reactions  of  Fe+ 
with  C2H 2,  C2H4,  and  C2H6  as  a  function  of  He  concentration.  Least 
squares  lines  forced  through  zero  are  shown. 


sextet  ground  state  of  the  Fe+  reactant.  In  Table  2,  we  list  both 
the  adiabatic  bond  dissociation  energy  (BDE)  (i.e.,  the  energy 


Table  2.  Calculated  BDEs  at  the  B3LYP/TZVP  Level,  as 
Well  as  the  Published  Experimentally  Determined  BDE  for 
Each  Association  Product” 


calculated  BDE  (eV) 


association  product 

adiabatic 

diabatic 

experimental  BDE  (eV) 

Fe(C2H2)+ 

2.14 

0.95 

1.S643'63 

Fe(C2H4)+ 

2.12 

1.00 

1.564 

Fe(C2H6)+ 

1.10 

0.34 

0.6641 

“The  adiabatic  BDE  assumes  spin  crossing  from  the  sextet  Fe+  to 
quartet  association  complex,  while  the  diabatic  BDE  assumes  no  spin 
crossing  (i.e.,  sextet  Fe+  and  sextet  association  complex). 


difference  between  the  quartet  association  complex  and  the 
sextet  Fe+  reactants)  and  diabatic  BDE  (i.e.,  between  the  sextet 
association  complex  and  the  sextet  Fe+  reactants).  For  all  three 


reactions,  previously  reported  experimental  BDEs  fall  about 
halfway  in  between  the  calculated  adiabatic  and  diabatic  BDEs. 

The  effective  rate  constants  reported  here  for  the  association 
of  Fe+  with  C2H2  and  C2H6  are  nearly  identical  over  the  entire 
temperature  range,  while  the  rate  constant  for  C2H4  is 
approximately  5  times  higher  and  declines  less  steeply  with 
temperature;  the  latter  is  probably  due  to  falloff  effects.  As 
previously  discussed  by  Weishar65  and  Bohme,44  the  relative 
difference  in  rate  constants  can  be  rationalized  in  terms  of  the 
complex  lifetime.  Fe(C2H2)+  and  Fe(C2H4)+  have  similar  bond 
dissociation  energies,43'63’64  but  because  Fe(C2H4)+  has  more 
degrees  of  freedom,  it  has  a  larger  density  of  states  at  a  given 
energy.  As  a  result,  calculating  the  unimolecular  dissociation 
rate  constants  of  the  complexes  via  orbiting  transition-state 
phase  space  theory  indicates  that  Fe(C2H4)+*  has  a  longer 
lifetime  than  Fe(C2H2)+*  and  thus  a  faster  rate  constant  for  the 
association.  More  accurate  calculations  need  to  know  the 
rigidity  factor,66  and  thus,  these  calculation  allow  only  for  a 
relative  comparison.  Fe(C2H6)+  has  the  most  complexity  but 
also  a  much  lower  BDE  (0.66  eV),41  leading  to  a  lifetime  similar 
to  that  for  Fe(C2H2)+  (at  the  dissociation  threshold).  The 
phase  space  theory  calculations  are  in  agreement  with  the 
observed  ordering  of  the  rate  constants.  The  relative  difference 
in  temperature  dependencies  is  in  part  related  to  the  reaction  of 
Fe+  with  ethylene  being  in  the  falloff  region;  the  high-pressure 
rate  constant  will  have  a  lower  temperature  dependence  than 
the  low-pressure  one. 

Accurate  modeling  of  the  observed  reactivity  in  this  case  is 
further  complicated  by  the  potential  effects  of  intersystem 
crossing  in  these  reactions.  So-called  two-state  reactivity,  where 
crossing  to  a  different  spin  potential  surface  allows  for  a  lower- 
energy  path  to  products,  has  been  observed  in  numerous 
reactions  involving  transition  metals,19’20  most  notably  with 
Fe+, 30,67-69  though  the  extent  of  the  role  that  it  plays  is  hard  to 
predict.  In  determining  the  BDE  of  Fe(C2H4)+,  Armentrout 
reports  enough  diabatic  (spin-conserving)  nature  to  the 
observed  thresholds  that  the  reported  BDE  is  adjusted  by  the 
4Fe+— sFe+  energy  splitting.64  Conversely,  we  have  recently 
reported  on  the  reaction  of  Fe+  with  N20,  finding  that  a 
significant  likelihood  of  curve  crossing  from  the  sextet  to 
quartet  surface  is  needed  to  model  the  observed  reaction  rate 
constants.70  Quantitative  understanding  of  the  interplay  of  spin 
states  in  reactions  of  this  type  is  still  in  its  infancy;  thus, 
attempts  to  model  this  system  in  greater  detail  would  be 
premature. 


Table  3.  Rate  Constants  (fc),  Product  Branching,  Collision  Rate  Constant  (fccou),  and  Reaction  Efficiency  (k/kCtt]])  for  the 
Bimolecular  Reactions  of  FeO+  with  Acetylene,  Ethylene,  and  Ethane  at  300  Ktt 


product  branching 


reaction 

this  work 

ut. 

Fe+0  +  C2H4 

Fe+  +  C2H20 

0.43 

o.sb 

->  FeCH2+  +  CO 

0.57 

o.sb 

FeO+  +  C2H4 

-►  Fe+  +  C2H40 

0.60 

0.82c 

-*■  FeCH,+  +  CH20 

0.26 

0.11c 

->  FeC2H2+  +  H20 

0.04 

0.04c 

->  FeCO+  +  CH4 

0.10 

0.03c 

FeO+  +  C2H6 

->  Fe+  +  C2H60 

0.12 

0.10, d  0.1 2e 

-►  FeC2H4+  +  H,0 

0.67 

0.70d  0.67e 

-►  FeH2+  +  C2H4 

0.21 

0.20/  0.1 2e 

k  (300  K)  cm3  s  1 

fccoii  cm3  s  1 

k/kcon  (300  K) 

T  dep 

9.8S  X  10~10 

1.06  X  10~9 

0.93 

~T° 

1.08  X  10-9 

1.08  X  10~9 

1.00 

~T° 

8.56  X  10  I0,  8.4  X  10~lod 

1.07  X  10~9 

0.80 

rp-0.47±0.15 

aErrors  in  the  rate  constants  are  estimated  to  be  ±25%  absolute  and  ±15%  relative  to  each  other.51  Errors  in  the  product  branching  are  estimated  to 
be  ±10%.  The  temperature  dependence  (T  dep)  is  also  given.  ^From  ref  13.  cFrom  ref  72.  ^From  ref  71.  eFrom  ref  45. 
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2.  FeO+  +  C2Hx.  The  rate  constants  and  product  branching 
observed  for  the  bimolecular  reactions  of  FeO+  with  acetylene, 
ethylene,  and  ethane  at  300  K  are  listed  in  Table  3,  as  are 
previously  published  values,  where  available. 13,45,7 1,72  The 
reaction  efficiencies,  defined  as  the  ratio  of  the  experimentally 
measured  rate  constant  to  the  collision  rate  constant  as 
calculated  by  the  parametrized  Su  and  Chesnavich  method,73 
are  shown  in  Figure  3  from  170  to  700  K.  Both  acetylene  and 


O 

2 


O 

c 

0 

O 


ifc 

LU 

c 

o 

•*-< 

o 

0 

0 


O' 


FeO+  +  C  H 

2  x 


100 


1000 


T(K) 


Figure  3.  Reaction  efficiencies,  defined  as  the  ratio  of  the 
experimentally  measured  rate  constant  to  the  calculated  collision 
rate  constant,  for  the  reaction  of  FeO+  with  acetylene,  ethylene,  and 
ethane  from  170  to  700  K. 


ethylene  react  at  the  collision  rate  over  the  entire  temperature 
range.  Ethane  reacts  at  the  collision  rate  at  170  K,  and  the 
efficiency  decreases  as  approximately  T~1/2.  Each  reaction  yields 
several  products,  which  are  discussed  below. 

The  behavior  of  these  reactions  is  entirely  consistent  with 
calculations  of  the  surfaces  that  have  been  published  elsewhere. 
Reaction  pathways  for  FeO+  with  acetylene47  and  ethane48  have 
been  published  by  Zhao  et  al.  employing  density  functional 
methods.  Sun  et  al.  further  calculated  the  energetics  of  DFT- 
optimized  stationary  points  on  the  reaction  surface  at  the 
CCSD(T)  level  for  the  reaction  with  ethane;  they  also 
introduced  a  potentially  competitive  direct  hydrogen  abstrac¬ 
tion  pathway  toward  the  oxygen-transfer  product  ethanol.42 
This  pathway  has  a  calculated  barrier  that  is  competitive  with 
the  concerted  pathway,  though  just  above  the  reactant  energy. 
To  the  best  of  our  knowledge,  no  surface  for  the  reaction  with 
ethylene  has  been  reported. 

For  the  acetylene  reaction,  the  calculated  highest-energy 
transition-state  barrier  for  the  lowest-energy  pathway  is  1.33  eV 
below  the  reactants,  while  the  products  lie  2.11  and  2.16  eV 
below.47  This  highly  exothermic  reaction  is  only  minimally 
hindered  kinetically  and  therefore  would  be  expected  to 
proceed  with  very  high  efficiency  and  minimal  variation  over 
this  temperature  range,  in  excellent  agreement  with  what  we 


observe.  While  no  surface  is  available  for  the  reaction  with 
ethylene,  a  similar  behavior  of  the  rate  constant  strongly 
suggests  that  the  largest  transition-state  barrier  lies  well  below 
the  initial  reactants  as  well.  For  the  reaction  of  FeO+  with 
ethane,  the  behavior  of  the  rate  constant  is  in  better  agreement 
with  the  surface  published  by  Zhao  et  al.48  In  this  case,  the 
largest  transition-state  barrier  is  0.55  eV  below  the  reactants, 
whereas  Sun  et  al.  found  the  largest  transition-state  barrier  on 
the  lowest-energy  pathway  to  be  only  0.069  eV  below  the 
reactants.42  In  both  cases,  traversing  the  lowest-energy 
transition-state  barrier  requires  crossing  from  the  sextet  surface 
to  the  quartet  surface,  potentially  an  example  of  two-state 
reactivity.19  Sun  et  al.  also  calculated  the  spin— orbit  coupling  at 
the  crossing  point  just  prior  to  the  0.069  eV  exothermic  barrier, 
suggesting  that  the  crossing  may  be  quite  facile.42  However,  a 
near-thermoneutral  barrier  such  as  this  is  inconsistent  with  a 
reaction  as  efficient  as  that  observed,  even  assuming  that  the 
crossing  ensues  with  no  kinetic  hindrance.  The  present 
measurements  are  therefore  more  in  line  with  the  pure  DFT 
calculations  than  the  DFT  optimizations  with  CCSD(T)  single¬ 
point  energies.  The  temperature  dependence  of  this  reaction  is 
explained  by  rotational  channel  switching.74  As  the  temperature 
rises,  both  the  internal  energy  distribution  and  the  /  distribution 
of  the  system  are  shifted  upward,  the  latter  due  to  the 
increasing  likelihood  of  large  impact  parameter  collisions 
surmounting  the  centrifugal  entrance  barrier.  Both  factors 
favor  returning  to  reactants  through  the  loose  entrance 
transition  state,  as  opposed  to  continuing  to  products  over 
the  tighter  transition  state,  the  rotational  barrier  for  which  rises 
more  steeply  with  /. 

The  product  branching  for  the  reactions  of  FeO+  with 
acetylene,  ethylene,  and  ethane  from  170  to  700  K  are  shown  in 
Figures  4—6,  respectively.  The  reaction  of  FeO+  with  ethane 
shows  strong  selectivity  (~67%)  for  production  of  FeC2H4+ 
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Figure  4.  Product  branching  for  the  reaction  of  FeO+  with  ethane 
from  170  to  700  K.  The  fits  shown  are  to  guide  the  eye. 
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Figure  5.  Product  branching  for  the  reaction  of  FeO+  with  ethylene 
from  170  to  700  K.  The  fits  shown  are  to  guide  the  eye. 
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Figure  6.  Product  branching  for  the  reaction  of  FeO+  with  acetylene 
from  170  to  700  K.  The  fits  shown  are  to  guide  the  eye. 


and  H20,  with  minor  production  of  FeH20+  and  C2H4  (~21%) 
and  very  little  production  of  Fe+  and  ethanol  (~12%).  The 
product  branching  remains  essentially  flat  over  the  entire 
temperature  range  studied.  This  result  agrees  with  the 
calculated  thermodynamics  for  the  reaction.  The  products  are 
well-separated  energetically,  lying  2.25,  1.65,  and  0.63  eV  below 
ground-state  reactants  as  calculated  by  Zhao  for  the  FeC2H4+, 


FeH20+,  and  Fe+  channels,  respectively.48  The  small  amount  of 
energy  added  with  temperature  compared  to  these  values  is 
consistent  with  the  small  dependence  of  the  branching  on 
temperature.  The  direct  hydrogen  abstraction  pathway  to  the 
Fe+  +  ethanol  channel  suggested  by  Sun  et  al.42  is  calculated  to 
have  a  rate-limiting  barrier  just  slightly  higher  than  that  for  the 
concerted  pathway  (0.078  to  —0.069  eV  relative  to  the 
reactants,  respectively),  and  the  curve  crossings  prior  to  these 
barriers  are  expected  to  be  comparable  (spin— orbit  coupling 
matrix  elements  of  315.5  and  306.6  cm-1,  respectively), 
indicative  of  this  pathway  likely  being  competitive.  As  discussed 
earlier,  however,  the  high  efficiency  of  this  reaction  is 
inconsistent  with  rate-limiting  barriers  near  thermoneutral, 
such  as  these.  Furthermore,  the  singular  temperature  depend¬ 
ence  of  the  rate  constant  and  especially  the  lack  of  temperature 
dependence  in  the  product  branching  strongly  support  a  single 
bottleneck  to  reactivity  and  thus  a  single  mechanism  over  this 
temperature  range.  While  the  direct  hydrogen  abstraction 
pathway  may  be  relevant  at  higher  energies,  near  thermal 
energies,  it  appears  to  have  no  impact. 

The  product  branchings  for  FeO+  with  ethylene  and 
acetylene  show  striking  similarities.  The  primary  product  ions 
formed  in  both  cases  are  FeCH2+  and  Fe+,  with  the  likely 
corresponding  neutrals  CFI20  and  C2H40  for  ethylene  and  CO 
and  C2H20  for  acetylene.  Two  additional  minor  products  are 
also  observed  for  the  reaction  with  ethylene,  FeCO+  +  CH4  and 
FeC2H2+  +  H20.  Both  of  these  minor  product  channels  remain 
small  over  our  temperature  range,  with  the  former  channel 
decreasing  slightly  and  the  latter  increasing  slightly.  The 
oxygen-transfer  channel  resulting  in  Fe+  formation  is  more 
prevalent  in  the  reaction  with  ethylene  than  that  with  acetylene 
over  all  temperatures  studied,  though  in  both  reactions,  this 
channel  decreased  substantially  with  temperature  while  the 
FeCH2+  channel  increased  significantly. 

Interpretation  of  the  product  branching  ratios  with  respect  to 
the  published  calculated  surfaces  is  difficult.  No  surface  has 
been  published  for  the  reaction  with  ethylene,  and  the 
calculated  surface  for  the  reaction  with  acetylene  is  extremely 
complicated,  with  numerous  intermediates  and  transition  states 
as  well  as  several  available  Fe+  channels  due  to  varied 
conformations  of  the  C2H20  neutral.47  Producing  Fe+  and 
neutral  CHCHO  is  calculated  to  be  1.19  eV  endothermic  and  is 
unlikely  to  be  involved  under  these  experimental  conditions. 
Changing  the  neutral  isomer  to  CHCOFI  makes  the  reaction 
0.57  eV  exothermic,  and  production  of  CH2CO  is  2.16  eV 
exothermic.  This  is  just  0.05  eV  below  the  FeCH2+  +  CO 
channel.  Further  complicating  matters  is  the  role  of  spin 
conservation  in  this  reaction.  As  previously  discussed,  two-state 
reactivity  is  prevalent  in  many  reactions  involving  transition- 
metal  cations.  The  calculated  surface  for  the  reaction  with 
acetylene  involves  numerous  transition-state  barriers  that  are 
considerably  lower  on  the  quartet  surface  as  opposed  to  that  on 
the  sextet  surface,  indicative  of  potential  two-state  reactivity. 
Furthermore,  the  ground  state  of  the  FeCH2+  product  is 
calculated  to  be  quartet,  while  the  ground  state  for  the  Fe+ 
channels  lie  on  the  sextet  surface,  strongly  suggesting  that 
intersystem  crossing  between  spin  states  has  significant 
influence  over  the  observed  product  branching. 

Because  a  primary  motive  for  this  study  is  furthering  our 
understanding  of  transition-metal  cations  in  real  world  catalysts, 
it  is  informative  to  consider  these  reactions  in  this  context.  As 
an  example,  we  consider  the  effectiveness  of  Fe+  as  an  oxygen- 
transfer  catalyst  for  these  three  hydrocarbons  in  combination 
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with  the  complementary  reaction  of  Fe+  +  N20,  which 
completes  the  catalytic  cycle  by  producing  FeO+  in  each 
reaction.  We  then  calculate  the  number  of  oxidized  products 
(alcohols)  produced  per  cycle,  accounting  for  products  that 
may  or  may  not  lead  to  completion  of  the  catalytic  cycle 
through  further  reaction  with  N20.  The  results  are  shown  in 
Figure  7.  An  ideal  catalyst  will  produce  one  oxygen-transfer 


Number  of  Cycles 

Figure  7.  The  relative  efficiency  of  the  oxygen-transfer  process  from 
163  to  600  K  for  the  reaction  of  FeO+  with  acetylene,  ethylene,  and 
ethane.  An  indication  of  an  ideal  catalyst  producing  one  oxygen- 
transfer  product  per  FeO+  per  cycle  is  shown  as  a  dashed  line.  Product 
channels  that  lead  to  Fe+,  either  initially  or  secondarily  with  further 
reaction  with  N20,  are  taken  to  resupply  FeO+  with  unit  efficiency. 


product  per  FeO+  per  cycle,  and  that  is  shown  as  the  1:1  line. 
The  bands  give  the  range  for  the  two  temperature  extremes, 
highlighting  the  variation  of  product  branching  for  the  reactions 
over  these  temperatures.  By  far,  the  most  efficient  reaction  is 
that  with  acetylene.  While  this  does  not  equal  the  ideal  case, 
because  FeCH2+  is  produced  as  a  side  product,  the  curve  does 
not  deviate  from  linear  with  increasing  cycles  because  FeCH2+ 
will  further  react  with  N20  to  resupply  Fe+,  and  thus,  the  cycle 
is  not  poisoned.  The  poorest  efficiency  is  seen  in  the  reaction 
with  ethane.  After  just  a  few  cycles,  the  number  of  alcohol 
molecules  produced  saturates,  indicating  that  no  more  are 
produced  due  to  the  poisoning  by  other  channels.  The  low 
turnover  number  for  this  reaction,  ~2.5,45  is  further 
compounded  by  the  poor  selectivity  toward  ethanol  produc¬ 
tion.  Oxygen  transfer  to  ethylene  is  an  interesting  case.  While 
the  reaction  is  more  selective  toward  oxygen  transfer  than  the 
reaction  with  acetylene,  the  seemingly  unimportant  minor 
channels  leading  to  Fe(CO)+  and  Fe(C2H2)+  provide  no  way 
back  to  Fe+  and  are  therefore  dead  channels.  This  limits  the 
turnover  number  for  this  reaction  to  ~4.8,  severely  hindering 
the  oxygen-transfer  efficiency. 

■  CONCLUSIONS 

We  have  studied  the  temperature  dependence  of  the  rate 
constants  and  product  branching  for  the  reactions  of  Fe+  and 
FeO+  with  acetylene,  ethylene,  and  ethane.  Fe+  is  observed  to 
react  with  the  three  hydrocarbons  only  by  association.  The 
relative  differences  in  reaction  rate  constants  are  due  to 
variations  in  the  complex  lifetime,  while  the  variations  in  the  T 
dependence  are  due  to  the  reactions  with  acetylene  and  ethane 


being  in  the  low-pressure  limit  for  ternary  reactions  while  the 
reaction  with  ethylene  is  in  the  falloff  regime.  The  reactions  of 
FeO+  with  the  three  hydrocarbons  are  very  efficient;  the 
reactions  with  acetylene  and  ethylene  take  place  at  the 
collisional  limit  over  the  entire  temperature  range  studied, 
and  the  reaction  with  ethane  is  collisional  at  170  K  and  falls  off 
at  T-0'5.  The  product  branching  for  the  reactions  of  FeO+  with 
acetylene  and  ethylene  both  showed  a  strong  T  dependence  in 
the  Fe+  and  the  FeCH2+  producing  channels,  with  the  latter 
being  favored  at  higher  temperatures  at  the  expense  of  the 
former,  while  the  product  branching  of  the  ethane  reaction 
remained  static  over  our  temperature  range. 

Each  of  the  studied  reactions  was  compared  to  published 
calculations  of  their  reactive  surfaces,  where  available.  Detailed 
kinetic  modeling  was  not  performed  due  to  the  poorly 
quantified  influence  of  two-state  reactivity  for  the  association 
reactions  and  extremely  complicated  nature  of  the  surfaces  for 
the  reactive  ones.  Qualitative  assessment  was,  however, 
possible.  The  high  efficiency  of  the  reactions  with  FeO+  was 
in  good  agreement  with  the  deeply  submerged  rate-limiting 
barriers  calculated  by  Zhang  et  al.  for  the  reactions  with 
acetylene  and  ethane,  with  the  latter  being  less  submerged,  in 
agreement  with  the  slightly  reduced  reactivity.  The  product 
branching  for  the  reaction  with  acetylene  showed  significant 
temperature  dependence  and  is  complicated  by  intersystem 
crossing  between  the  quartet  and  sextet  potential  energy 
surfaces.  The  product  branching  for  the  ethane  reaction  is  flat 
over  our  temperature  range,  consistent  with  the  well-spaced, 
highly  exothermic  product  channels  calculated  by  Zhang  et  al. 
The  static  product  branching,  combined  with  the  high 
efficiency  of  the  reaction  as  well  as  a  single  temperature 
dependence,  suggests  that  the  alternate  reaction  pathway 
calculated  by  Sun  et  al.  has  little  impact  at  thermal  energies. 

The  comparison  of  the  oxygen-transfer  product  channel  for 
the  reactions  of  FeO+  with  the  three  hydrocarbons  highlights 
the  importance  of  understanding  the  variation  in  product 
branching  for  these  types  of  reactions.  The  reactions  with  both 
acetylene  and  ethylene  exhibit  a  strong  temperature  depend¬ 
ence  to  their  respective  product  branching  and  hence  the 
selectivity  of  the  reaction.  More  importantly,  however,  are 
channels  that  do  not  allow  for  a  recycling  of  the  catalytic  ion. 
While  the  reaction  of  FeO+  with  ethylene  shows  more 
selectivity  toward  oxygen  transfer  than  acetylene  at  all 
temperatures  studied,  the  reaction  with  acetylene  is  significantly 
more  productive  as  none  of  the  products  poison  the  catalytic 
cycle. 
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